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Abstract 

Standard active-sterile neutrino oscillations do not provide a satisfactory description of the 
LSND evidence for neutrino oscillations together with the constraints from MiniBooNE and other 
null-result short-baseline oscillation experiments. However, if the mass or the mixing of the sterile 
neutrino depends in an exotic way on its energy all data become consistent. I explore the phe- 
nomenological consequences of the assumption that either the mass or the mixing scales with the 
neutrino energy as l/i?^ (r > 0). Since the neutrino energy in LSND is about 40 MeV, whereas 
MiniBooNE operates at around 1 GeV, oscillations get suppressed in MiniBooNE and the two 
results become fully compatible for r ^ 0.2. Furthermore, also the global fit of all relevant data 
improves significantly by exploring the different energy regimes of the various experiments. The 
best fit decreases by 12.7 (14.1) units with respect to standard sterile neutrino oscillations if 
the mass (mixing) scales with energy. 
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I. INTRODUCTION 



Reconciling the LSND evidence pL|] for t'u — oscillations with the global neutrino data 
reporting evidence 

mm 

and bounds 

mm 

, [lol, [ll| on oscillations remains a long- 
standing problem for neutrino phenomenology. Recently the MiniBooNE experiment 12 
added more information to this problem. This experiment searches for —>■ v^. appearance 
with a very similar LjEy range as LSND. No evidence for oscillations is found and the results 

The 



12 



are inconsistent with a two-neutrino oscillation interpretation of LSND at 98% CL 
standard "solution" to the LSND problem is to introduce one or more sterile neutrinos at 
the eV scale [l3|]. However, it turns out that such sterile neutrino schemes do not provide a 
satisfactory description of all data in terms of neutrino oscillations, see Ref. ij] for a recent 



analysis including MiniBooNE data; pre-MiniBooNE analyses can be found, e.g., in Refs. 
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Apart from sterile neutrino oscillations, various more exotic explanations of 

CPT violation 



the LSND signal have been proposed, for example, neutrino decay 
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, violation of Lorentz symmetry 
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CPT-violating quantum decoherence 
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varying neutrinos 25l| . or shortcuts of sterile neutrinos in extra dimensions 26l |. See Refs. 271 . 
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mass- 



for two recent reviews. 
In view of the difficulties to describe all data with "standard" sterile neutrino oscillations 
I assume in this note that the sterile neutrino is a more exotic particle than just a neutrino 
without weak interactions. Being a singlet under the Standard Model gauge group it seems 
possible that the sterile neutrino is a messenger from a hidden sector with some weired 
properties. I will assume in the following that the mass or the mixing of the fourth neutrino 
depends in a non-standard way on the neutrino energy. The motivation is that the neutrino 
energy in LSND is around 40 MeV, whereas MiniBooNE and the CDHS disappearance 
experiment operate around 1 GeV, and hence changing the energy dependence of oscillations 
will have some impact on the fit. 

To be specific, I am going to assume that either the mass m4 or the mixing [7^4 [a = e, /x) 
of the fourth neutrino state depends on the neutrino energy Ei^ like 

Erei 



MED: 



EDM: 



m^{Ei, 



ml 



E,. 



\U^4\\E,) = \U, 



a4 



E, 



ref 



E. 



(1) 



where m^ and lfa4. are the mass and the mixing at a reference energy Ej-^f, and for the 
exponent r I will consider values in the interval < r < 1. An energy dependent mass as in 
Eq. ([1]) will modify the energy dependence of oscillations from the standard l/E^, dependence 
to 1/El'^^'. Therefore, I refer to this effect as modified energy dependence (MED) oscillations, 
whereas the second case is denoted by energy dependent mixing (EDM). Most likely rather 
exotic physics will be necessary to obtain such a behaviour, some speculative remarks on 
possible origins are given in Sec. El Here I make the assumptions of Eq. ([T]) without 
specifying any underlying model, instead I will explore the phenomenological consequences 
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of MED and EDM oscillations for short-baseline neutrino data. I will show that in both cases 
{i) LSND and MiniBooNE become compatible and (ii) the global fit improves significantly. 

It is important to note that Eq. ([1]) involves only the new fourth mass state, whereas 
masses and mixing of the three active Standard Model neutrinos are assumed to be constant. 
Therefore the successful and very robust description of solar, atmospheric, and long-baseline 
reactor and accelerator experiments 

mm 

by three-flavour active neutrino oscillations 
is not altered significantly, since the mixing of the forth mass state with active neutrinos is 
small. 

The outline of this paper is as follows. In Sec. [ITl I discuss in some detail the framework 
of MED and EDM oscillations and give a qualitative discussion of the expected behaviour of 
the combined analysis of the relevant short-baseline oscillation data. In Sec. IIIII the results 
of the numerical analysis are presented. The global fit includes the appearance experiments 
LSND, MiniBooNE, KARMEN and NOMAD, as well as various disappearance experiments. 
In Sec.|IV]I comment briefly on phenomenological consequences of MED /EDM oscillations in 
future oscillation experiments, astrophysics and cosmology. Sec.|V]presents some speculative 
thoughts on models leading to energy dependent masses and mixing for sterile neutrinos, 
and I summarize in Sec. IVII 



II. THE MED AND EDM OSCILLATION FRAMEWORKS 



Before discussing qualitatively the phenomenological consequences of MED and EDM 
oscillations according to Eq. ([1]) let me briefly remind the reader about the description 
of short-baseline (SBL) neutrino oscillation data in the case of standard four-neutrino os- 
cillations. In the so-called (3+1) schemes there is a hierarchy among the mass-squared 
differences: 

Am2i < |Am3i| < Aml^ . (2) 

Under the assumption that Amgi and Amg^ can be neglected SBL oscillations are described 
by two-flavour oscillation probabilities with an effective mixing angle depending on the 



elements of the lepton mixing matrix |f/e4p and \Uf^i\ 
the effective mixing angle is given by 



For Ufj, — > z/e appearance experiments 



sin^ 26 



fie 



Ml 



(3) 



whereas for a z/q, disappearance experiment we have 

sin2 2^„„ = 4|f/„4|'(l-|t/a4|)' 



(4) 



see, e.g. 



Ref. 



29 



is a product of |[/e4p and |f/^4 



The fact that the amplitude responsible for the LSND appearance 
^, whereas I'e and i^^ disappearance experiments constrain 



these elements separately leads to the well-known tension between LSND and disappearance 



experiments in the (3+1) oscillation schemes, see e.g., Refs. [29, 



30 



31 



33 



34 
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Assuming now an energy dependent mass for 1/4 as in Eq. ([T]), it turns out that for the 
range of parameters and energies relevant for our discussion it is always possible to take 

much heavier than the three standard neutrinos, mi^2,3 ^4, such that the usual SBL 
approximation Eq. ([2]) remains always valid, and the energy scaling of Eq. ([T]) applies also 
for the mass-squared difference Am\i{Ey) = m\{Ey) — m\ ^ m1{E^). Then the relevant 
oscillation phase 0osc gets a different energy dependence than in the standard case: 

Hence the standard l/E^, dependence gets altered to 1/El'^^'. The most relevant consequence 
of the MED with r > follows from Eq. ([5]): An experiment is sensitive to oscillations if 

0OSC ^ vr/2, or 

L 2n f E,.y 

(6) 



E^ Amj-^ V-^rof, 

Hence, for experiments with Ei, > i^ref the allowed region will be shifted to larger values of 
Am|^ as compared to the standard oscillation case (r = 0), whereas the allowed region for 
experiments with E^ < i^ref will be shifted to smaller Arhli, in order to compensate for the 
factor {EjE,efY in Eq. (E]). 

Using instead of the MED now the energy dependence of the mixing matrix elements 
from Eq. ([T]), one obtains for the SBL appearance and disappearance amplitudes given in 
Eqs. (I3D and dH): 

/E \ 

sin2 2^^, = 4|[/e4n?7^4|'cx '''^ 



sin2 2^„„ = 4|f/„4p(l - |f/a4n ~ 4|f/,4p oc 



(7) 



This introduces only a mild distortion of the oscillation pattern from the standard oscillatory 
behaviour with l/Eu. The main effect of EDM is that the sensitivity of experiments with 
El, > Ej-cf to the effective mixing angle gets weaker. Let us note that the EDM scaling of 
Eq. ([1]) cannot hold for arbitrarily low energies, simply because of unitarity of U. The low 
energy limit of the EDM scaling should find an explanation in some theory for this effect. 
Here I assume that for the parameter range relevant for the SBL analysis the power-law 
scaling of Eq. ([T]) remains valid. 

Note that in both cases, MED and EDM, for fixed r the choice of the reference energy 
Eref is arbitrary. From Eq. ([5]) it is clear that choosing a different reference energy E^et leads 
just to a rescaling of m4 such that the combination Aml^E^^^ remains constant. Similar, 
changing Ej-et in case of the EDM leads just to a rescaling of the |f/Q,4p. Hence, the (3+1) 
MED and EDM models have one phenomenological parameter in addition to the (3+1) 
standard oscillation model: 

MED: \UeA\MUf,4\^ Arhl^, r , 
EDM: \Ue4\^,\U^,\^, Aml„r. 
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Experiment 


Ref. 


Channel 


Data 




Bugey 


[9J 


Uf. - 




60 


4 MeV 


Chooz 


[10] 


Vf. - 


-> 


1 


4 MeV 


Palo Verde 


M 






1 


4 MeV 


LSND 


in 


- 


Ve 


11 


40 MeV 


KARMEN 


[6] 


- 


^ Ve 


9 


40 MeV 


MiniBooNE 


[12j 


- 


^ Ve 


8 


700 MeV 


CDHS 


m 


- 


Vy. 


15 


1 GeV 


NOMAD 


m 


- 


^ Ve 


1 


50 GeV 



TABLE I: Experiments used in the SBL analysis. The oscillation channel, and the number of data points 
used in the fit ("Data") are given. The last column shows the approximate mean neutrino energy for each 
experiment. 

The main effects of MED and EDM oscillations can be summarized in the following 
way: Consider the allowed regions of the various experiments for standard oscillations in 
the plane of sin^ and Am|^. Introducing now MED (EDM) leads to a relative shift 
of the regions of experiments at different energies along the Am^j^ (sin^ axis. The 
relevant SBL experiments are listed in Tab. [H ordered according to their mean neutrino 
energy. For convenience I will choose E'ref = 40 MeV, corresponding roughly to the mean 
neutrino energy in LSND. Since in MiniBooNE the neutrino energy is higher, for MED 
(EDM) oscillations the sensitivity is shifted to larger values of ISm\^ (sin^2^^e) for r > 0. 
In the MED framework MiniBooNE operates actually at a too small value of L/ in order 
to test LSND. As I will show in the following, in both cases the two experiments become 
compatible for r > 0. Furthermore, it turns out that also the global fit including all the 
experiments listed in the table improves significantly due to the different energy regimes. 



III. GLOBAL ANALYSIS OF SBL DATA FOR MED AND EDM OSCILLATIONS 
A. Description of the data used in the fit 

Before presenting the results of the analysis let us briefly discuss the data used in the 
fit, as summarized in Tab. [B For the re-analysis of LSND I fit the observed transition 
probability (total rate) plus 11 data points of the L/E spectrum with free normalisation, 
both derived from the decay-at-rest data [l|. For KARMEN the data observed in 9 bins of 
prompt energy as well as the expected background |6|_is used in the fit. Further details of 
the LSND and KARMEN analyses are given in Ref. 
using the information provided in Ref. |7|. 
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For NOMAD I fit the total rate 



The MiniBooNE analysis is based on the information provided at the web-page 



36 



derived from the actual Monte Carlo simulation performed by the collaboration. Using these 
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data the "official" MiniBooNE analysis 12| can be reproduced with very good accuracy. 
The averaging of the transition probability is performed with the proper reconstruction 
efficiencies, and detailed information on error correlations and backgrounds is available for 
the analysis. MiniBooNE data are consistent with zero (no excess) above 475 MeV, 
whereas below this energy a 3. 60" excess of 96 ± 17 ± 20 events is observed. Whether this 
excess comes indeed from i/^ — > transitions or has some other origin is under investigation. 
As discussed in Ref. 



12 



standard two-neutrino oscillations cannot account for the event 
excess at low energies. Following the strategy of the MiniBooNE collaboration the analysis 
is restricted to the energy range from 475 MeV to 3 GeV. In Sec. IIII CI I will comment on 
the possibility to obtain the low energy event excess in case of MED oscillations. 
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reac- 



I include the disappearance experiments Bugey [9j, Chooz [10|, and Palo Verde 
tor Vf, disappearance), as well as the CDHS disappearance experiment, see Ref. 34 1 for 
technical details. In addition to the data listed in Tab. [T] atmospheric neutrino data give an 
important constraint on the mixing of z/^ with the heavy mass state, i.e., on |?7^ 



Ml 



32, 



35 



I use the updated analysis described in detail in Ref. [1J|. The basic assumption in this 
analysis is that Am^^ is "infinite" for atmospheric neutrinos according to Eq. ([2]). Since I 
assume this to be true also in the scenarios under consideration one can directly apply the 
bound on |t/^4p from the standard oscillation analysis. In case of EDM one should perform 
a re-analysis of atmospheric data taking into account the energy dependence of |?7^4p for the 
various data samples spanning five decades in neutrino energy. Such an analysis is beyond 
the scope of the present work and I assume a scaling of |t/^4p corresponding to an average 
energy of 1 GeV. Adding one data point for the bound from atmospheric neutrinos to the 
data given in Tab. [T]I obtain A^tot = 107 data points in the global analysis. 



B. Results of the global analysis 

Let us now discuss the results of the numerical analysis within the frameworks of MED 
and EDM oscillations. Figs. [T] and [2] show the allowed regions in the plane of t^m^^ and 
sin^ IQ^e for various data sets for standard oscillations (r = 0) compared to the exotic en- 
ergy dependence models for some values of r > 0. First, note that the neutrino energy 
in LSND and KARMEN is the same, and therefore the consistency of LSND and KAR- 
MEN 33] is not affected by introducing a non-zero r. Furthermore, since we choose a 
reference energy E'ref close to the mean energy in these experiments the allowed region from 
LSND+KARMEN-I-NOMAD does not change by increasing r} Second, in agreement with 
the argument given in Sec. [ITlone observes from the figures that the bound from MiniBooNE 
moves to higher values of lS.m\-^ for MED and to higher values of sin^ 26'^e for EDM if r in- 
creases, as a consequence of the higher neutrino energy in MiniBooNE. I find that in both 
cases for r ^ 0.2 LSND and MiniBooNE are fully consistent. Third, Figs. [1] and [2] show the 



^ NOMAD, with an energy of about 50 GeV, contributes very little to these regions. 
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A sterile neutrino with an energy dependent mass 




-4 -3 -2 -1 -4 -3 -2 -1 -4 -3 -2 -1 

10 10 10 10 10 10 10 10 10 10 10 10 

sm^2e 



FIG. 1: (3+1) oscillations with a modified energy dependence (MED): Allowed regions at 99% CL (2 d.o.f.) 
for LSND+KARMEN+NOMAD, MiniBooNE, and the disappearance experiments for standard oscillations 
r — (left) and MED oscillations with r = 0.3 (middle) and r = 0.6 (right). The star in the middle panel 
corresponds to the global best fit point. 



A sterile neutrino with energy dependent mixing 




-4 -3 -2 -1 -4 -3 -2 -1 -4 -3 -2 -1 

10 10 10 10 10 10 10 10 10 10 10 10 

single @ 40 MeV 



FIG. 2: (3+1) oscillations with energy dependent mixing (EDM): Allowed regions at 99% CL (2 d.o.f.) 
for LSND+KARMEN+NOMAD, MiniBooNE, and the disappearance experiments for standard oscillations 
r = (left) and EDM oscillations with r = 0.3 (middle) and r — 0.74 (right). The star in the right panel 
corresponds to the global best fit point. 

bound on sin^ 29^f, from the disappearance experiments Bugey, Chooz, Palo Verde, CDHS, 
and atmospheric neutrino data, where for a given sin^ 2^^e I minimize the with respect to 
|[/e4p and |t/;i4p under the constraint sin^ 2^^e = 4|?7e4p|t/^4p. As visible in the left panels, 
in the standard oscillation case there is severe tension between these data and the LSND 
evidence, and at 99% CL there is basically no overlap of the allowed regions. However, the 
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situation clearly improves for MED and EDM oscillations, and for r > the allowed regions 
overlap. 

For MED oscillations (Fig. [1]) this can be understood in the following way. The pro- 
nounced wiggles in the disappearance bound visible in the left panel around Aml^ ~ 1 eV^ 
come from the Bugey reactor experiment. Since Ei, for Bugey is smaller than the reference 
energy of 40 MeV these features move to lower values of Arhli if r is increased. On the other 
hand, the constraint on |?7^4p from CDHS is shifted to higher values of Am|j^ and only the 
weaker constraint from atmospheric data remains. Both trends work together in moving the 
disappearance bound towards the LSND region, as visible in the middle panel for r = 0.3. 
If r is further increased the Bugey pattern moves to even smaller values of Am||, and the 
bound at the relevant region around 1 eV^ is given by the constraints from Chooz on |t/e4p 
and atmospheric neutrinos on |f/^4p (see right panel). 

In the case of EDM oscillations (Fig. [2]) there are two opposite trends. Since for Bugey 
El, < i^ref the constraint on |f/e4p becomes stronger with increasing r, whereas for CDHS 
and atmospheric data E,j > E^^f and the bound on |f/^4p gets weaker. The upper limit on 
sin^26'^e emerges from the product of these two constraints, see Eq. ([2]), and therefore, it 
scales according to 



^ref 



Bugey (-E'i/)cDHS,atm 



0.4" . (9) 



Hence, the net-effect is a shift of the bound to larger values of sin 26'^e and a significant 
overlap with the LSND region emerges. 

Fig. [3] shows the for appearance data only (left) and for the global data (right) as a 
function of the energy exponent r. The best fit points have the following values: 



MED: x^pp,„.in = 19.4/(29 - 3) d.o.f. , xiob.min = 89.0/(107 - 4) d.o.f. 
EDM: xL,min = 19.6/(29 - 3) d.o.f. , xjob.min = 87.6/(107 - 4) d.o.f. 

and occur at the parameter values 



(10) 



MED: 




= 0.15, 




= 0.21, Am^i 


= 1.0 eV^ r = 


: 0.3, 


EDM: 


\ije4\ 


= 0.06, 




= 0.53, Aml-^ 


= 0.92 eV^ r : 


= 0.74 



(11) 

As mentioned above, the parameters with tilde correspond always to a reference energy of 
40 MeV. In agreement with the discussion above one observes that for EDM {Ueil is rather 
small to respect the stronger bound from Bugey, whereas |f/^4| gets relatively large due 
to the relaxed bound from CDHS and atmospheric data. Although the best fit point for 
EDM occurs at the relatively large value of r = 0.74, one can see from Fig. [3] that fits of 
comparable quality are obtained already for r ^ 0.4. 

The MED and EDM fits improves significantly with respect to the standard oscillation 
case r = 0: 



MED: Axl,{r = 0) = 7.7 , Ax^^^r = 0) = 12.7 . 
EDM: Ax^pp(r = 0) = 7.5 , Axl^r = 0) = 14.1 . 



(12) 



8 



35 



30 



25 



20 



15 



























































_ 










/appeal iiucc uai 


a 










_ 






































_ 






































_ 

























EDM 








_ 
























I 


vlEI 


) 








_ 

_ 






































_ 


















































































































































































































































































































































































) 


0.2 


0.4 


0.6 




0.8 



85 



105 



100 



95 































































lata 














_ 
_ 
















c 






















































_ 








































_ 















































































_ 








































_ 

_ 
















































































_ 














































\\ 














































































\ 








































\ 







































































































































































































































0.2 



0.4 0.6 
r 



0.8 



FIG. 3: The ^ for SBL appearance data (left) and global data (right) for MED and EDM oscillations as 
a function of the exponent r, minimized with respect to the other parameters. The number of data points 
is 29 for the appearance experiments and 107 for the global analysis. 



For comparison, the extension of the standard (3+1) oscillation scheme to a (3+2) scheme by 
the addition of a second sterile neutrino leads to an improvement of Xmin (3+1) ~ Xmin (3+2) ~ 
6.1 ij]. Taking into account that for (3+2) oscillations 4 additional parameters are intro- 
duced in the fit instead of only one as in the cases of MED or EDM, one concludes that the 
latter provide a much more significant improvement of the fit. 

The shape of the curves in Fig. [3] can be understood from the behaviour of the allowed 
regions shown in Figs. [T] and [21 The best fit for appearance data is reached once the 
MiniBooNE exclusion curve is moved out of the LSND region, and no further improvement 
can be obtained by further increasing r. In the case of EDM the fit gets worse again due to 
the energy distortion introduced for large r by scaling the mixing with 1/-E^. This is also 
the reason for the change in the LSND region visible in the right panel of Fig. [2l For MED, 
the global reaches a minimum when the wiggles from the Bugey bound cover the LSND 
region around l^^m\^^ ~ 1 eV^. If r is further increased these wiggles are moved out again of 
the LSND region and the fit gets slightly worse again. For r > 0.5 a plateau is reached, since 
then the disappearance bound at Am4i ~ 1 eV^ comes mainly from Chooz and atmospheric 
data, which are independent of lS.vn\^ and hence also independent of r. In the case of EDM, 
the global improves until relatively large values of r as a consequence of Eq. (jH]). At some 
point again the fit gets worse due to the anomalous energy dependence of the probability. 

A powerful tool to evaluate the compatibility of different data sets is the so-called para- 
meter goodness-of-fit (PG) criterion discussed in Ref. j38|. It is based on the function 



ot,min 



2 

i,min ' 



(13) 



f j^in is the minimum of the 



where Xtot,min ^^e minimum of all data sets combined and x 
data set i. This x^ function measures the "price" one has to pay by the combination of the 
data sets compared to fitting them independently. It should be evaluated for the number of 
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Standard oscillations MED oscillations EDM oscillations 
Data sets Xpc/d-o.f. PG xlc/d-O-f- PG xlc/d-o-f- PG 

LSND vs NEV 24.9/2 4 x lO'^ 14.0/3 0.3% 11.9/3 0.8% 

LSND vs NEV-APP vs DIS 25.3/4 4 x IQ-^ 14.3/6 3% 12.3/6 5% 
LKN vs MiniBooNE vs DIS 20.1/4 5 x 10"^ 8.9/6 18% 6.7/6 35% 

TABLE II: Consistency tests of various data sub-sets for standard (3+1) oscillations, MED and EDM 
oscillations. The data sets are SBL data showing no evidence for oscillations (NEV), no-evidence appearance 
data (NEV-APP = MiniBooNE + KARMEN + NOMAD), SBL disappearance data (DIS), and LSND + 
KARMEN + NOMAD (LKN). 1 give Xpg according to Eq. ^ and the corresponding probability ("PG"). 
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d.o.f. corresponding to the number of parameters in common to the data sets, see Ref. 
for a precise definition. 

The results of such a PG analysis are displayed in Tab. |IIl First, the compatibility of 
LSND and all the remaining no-evidence SBL data is tested, and the PG is compared within 
the standard, the MED, and the EDM oscillation frameworks. The consistency improves 
drastically from 4 x 10"^ to 3 x 10"^ (MED) or 8 x 10"^ (EDM). The probability value in 
the MED case corresponds to a tension of slightly less than 3a. Second, as an alternative 
test I check the compatibility of the three data sets LSND, no-evidence appearance data, 
and disappearance data. Similar a huge improvement of the consistency from 4 x 10~^ for 
standard oscillations to 3% (MED) or 5% (EDM) is found. 

Although both of these two tests clearly show an improvement of the fit with respect to 
standard oscillations, the low probabilities still indicate that the global fit is not perfect and 
some tension remains in the data. The reason for this is the tension between KARMEN and 
LSND, which is the same as in the standard oscillation case, since these experiments have 
the same energy. The allowed regions in the space of oscillation parameters of these two 
experiments have clearly an overlap, and a careful combined analysis came to the conclusion 



that they are consistent [37[. Nevertheless, there remains a tension between them which is 
detected by the rather sensitive PG test. In the last row of Tab. [IT] I consider the case when 
LSND and KARMEN are combined to one single data set (which includes also NOMAD), 
and test the consistency of this set against MiniBooNE and disappearance data. This 
corresponds to the data sets shown in Figs. [1] and [21 and since LSND and KARMEN are 
included in the same data set the tension between them does not show up in the PG. In 
this case the PG shows a perfect consistency of all data with probablities of 18% and 35% 
for MED and EDM, respectively, whereas the probablitiy of standard oscillations remains 
unacceptably low. 



10 



C. The low energy excess in MiniBooNE 



Before concluding this section I comment briefly on the event excess observed in Mini- 
BooNE in the energy region below 475 MeV. As discussed in Ref . 1^ , standard two-flavour 
oscillations cannot account for the sharp rise at low energy. However, since in the MED 
scenario the energy dependence of oscillations is modified according to Eq. 1^ one may 
expect that an explanation of the excess becomes possible. Indeed I find that for values 
of the exponent r ^ 1 the rise of the oscillation probability becomes steep enough and a 
perfect fit to the full spectrum including the excess between 300 and 475 MeV becomes 
possible. For such large values of r a closed allowed region appears in the plane of Am|^ and 
sin^ 29 for MiniBooNE data (not only a bound). However, because of the large r value this 
allowed region appears at Am|^ values above the LSND region and the KARMEN bound. 
Hence, although a modified energy dependence like considered here allows in principle for 
an explanation of the low energy event excess, this solution is not compatible with LSND 
and the KARMEN bound, and in the global analysis the excess cannot be fitted. For this 
reason I used in the analysis only the MiniBooNE data above 475 MeV 13], and rely on an 
alternative explanation of the low energy event excess. 



IV. IMPLICATIONS FOR FUTURE EXPERIMENTS, ASTROPHYSICS AND 
COSMOLOGY 

In this section I briefely comment on other phenomenological implications of the 
MED/EDM schemes. In general the scenarios considered here are difficult to test at fu- 
ture neutrino oscillation experiments. No appearance signal is expected for MiniBooNE 
anti-neutrino data, which currently are being accumulated, since the energy dependence is 
assumed to affect anti-neutrinos in the same way as neutrinos. In order to test the LSND 
signal for a MED with r = 0.3 at the given MiniBooNE baseline one would need to run at 
an energy of 

<fn?BooNE - ^ E,s.n - 40 MeV - 360 MeV , (14) 

V -^LSND / V 30 m / 

which seems not practicable because of low cross sections and large backgrounds. A similar 
signal as in standard (3+1) or (3+2) oscillations is expected also in the MED scenario for 



future reactor experiments, see Ref. [39|. A promising place to look for effects of MED 
oscillations could be the 2 km detector proposed for the T2K experiment jio!] . With a mean 
neutrino energy of 0.7 GeV this detector is slightly too far from the neutrino source to cover 
the oscillation maximum in case of standard oscillations with Am^ ~ 1 eV^. However, with 
the MED best fit point from Eq. ( fTTI) the oscillation phase according to Eq. ([5]) turns out 
to be close to 7r/2 at L = 2 km. 

A rather model independent test of MED or EDM explanations of LSND would be an 



experiment operating at the same energy as LSND, such as proposed in Ref. 
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Further- 
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more, to test the EDM scenario one would like to perform experiments at energies as small 
as possible. In particular, this model predicts relatively large effects for i/^ disappearance 
for experiments at energies smaller than 1 GeV, since the eV-scale mass state has a rather 
large mixing with z/^ at 40 MeV, see Eq. ([TT]) . Let us note that unitarity requires \Ua4,\ < 1- 
Therefore, the power law energy dependence of EDM cannot hold down to arbitrarily low 
energies. At the best fit values given in Eq. ffTTl) . one finds \Ue4\ ~ 1 at -Ej^ ~ 0.02 MeV 
and \Ufj_4^\ ~ 1 at -El, ~ 7 MeV.^ In order to use the EDM framework for very low energies 
one would have to specify the energy dependence of the neutrino mass matrix, and obtain 
the mixing angles via the diagonalisation, such that unitarity is always guaranteed. This is 
especially relevant, for example, to obtain predictions for neutrino mass experiments from 
Tritium beta-decay, which has an end point energy of 18.6 keV. 

As in case of standard sterile neutrino mixing, also in the MED/EDM framework the ster- 
ile neutrinos have implications for cosmology and astrophysics, see Ref. ^3] and references 
therein. In general the effects will be very similar to the standard case with effective masses 
and mixing evaluated according to the relevant neutrino energy. For example, in a supernova 
and in Big Bang nucleosynethesis (BBN) the neutrino energy is close to (or slightly below) 
the reference energy 40 MeV used above, and therefore the mixing parameters shown in 
Eq. (fTT!) roughly apply in these environments. This implies that — as in the standard (3-1-1) 
case — the sterile neutrino will be brought into thermal equilibrium via oscillations prior to 



BBN 



30, 
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Cosmology provides a bound on the sum of the neutrino masses in the sub-eV range 
mainly from the power spectrum at large scales combined with precise data on the cosmic 
microwave background. In general this bound implies a challenge for sterile neutrino schemes 
relevant for LSND. The conflict becomes particularly severe for the MED framework, since 
here the neutrino mass increases with decreasing neutrino energy, which implies large masses 
for cosmological relevant neutrinos. Let us note, however, that depending on the particular 
model realisation of MED one can expect that the power law scaling of Eq. ([1]) does not 
continue down to arbitrarily low energies. Here I assume only that it holds in the energy 
interval relevant for SBL experiments, i.e., above about 1 MeV, and it might be altered at 
lower energies. 



V. SPECULATIONS ON MODEL REALISATIONS OF MED/EDM 

Before concluding I give here some speculative thoughts on possible reasons for a power 
law scaling of sterile neutrino masses or mixing. Without doing any detailed model building 
I just mention a few possibilities where such a behaviour might occur. 

Indeed, energy dependent neutrino masses and mixing are a very familiar phenomenon in 



^ Note that for the SBL analysis only \Ue4\ is needed at few MeV energies, whereas \Uf^4\ is evaluated only 
for E,, > 40 MeV. 
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the framework of the standard matter effect 45] . Since the effective matter potential depends 



on the neutrino energy the mass eigenstates and mixing angles in matter depend on the 
energy. An analogous mechanism would be at work for the sterile neutrino if it interacts with 
some un-known background field or has some special interactions with standard matter. Such 



a possibility has been noted in Ref. [27| and explored recently in Ref. [46|. The interaction 
postulated for the sterile neutrino should be several orders of magnitude stronger than usual 
weak interactions in order to be relevant at the short baselines in LSND or MiniBooNE. In 
these models neutrinos and anti-neutrinos interact differently and therefore they will have a 
different energy dependence. This may change the null-prediction for the MiniBooNE search 
with anti-neutrinos mentioned in the previous section. 
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An energy dependence similar to a matter potential occurs also in the model of Ref. 
where sterile neutrinos are allowed to take shortcuts through particularly shaped extra 
dimensions. Effectively this leads to a modification of the dispersion relation of the sterile 
neutrinos which introduces a non-standard energy dependence on active-sterile oscillations. 
In general also a violation of the Lorentz symmetry such as considered for example in Ref. 23 
leads to oscillations with an energy dependence different from the standard one. 

Another motivation for energy dependent neutrino masses and/or mixing could be the 
idea of "unparticle" physics [47|. One assumes the existence of a scale invariant sector with a 
non-trivial infrared fixed point coupled to the Standard Model through non-renormalizable 
operators. Such operators may have large anomalous dimensions and hence introduce power- 
law running of coupling constants. For example, suppose a fermionic unparticle operator On 
with mass-dimension d^, with 3/2 < (In < 5/2, which has the quantum numbers of a right- 
handed neutrino (and hence is a gauge singlet). Then one can write a "mass term" rriaUaON, 
where z/„ can be either an active left-handed neutrino, or a "standard" sterile neutrino. The 
dimension of rria is 5/2 — with the anomalous dimension 3/2 — d^- Assuming that 
the coupling of On with is scale invariant implies that the effective "masses" at two 
energy scales Ai and A2 are related by mQ,(Ai) = (A2/Ai)'^'^~^/^mQ,(A2). This resembles 
the scaling of Eq. ([T]), relating the phenomenological parameter r in the above analysis to 
the anomalous dimension of the unparticle operator. The physical neutrino masses have to 
be found as poles in the corresponding propagator. These arguments provide a hint that 
the unparticle framework might lead to the exoting energy dependence of Eq. ([1]); whether 
it is indeed possible to obtain a valid model for MED and/or EDM active-sterile neutrino 
oscillations using unparticles needs further investigation, which is beyond the scope of this 
work. 

Via the so-called AdS / CFT correspondence effects from a conformal sector such as men- 
tioned above might actually have an interpretation also in theories with extra spacetime 

If the neutrino 



dimensions. In such models couplings can exhibit power law running 48 
mass is generated through a mechanism involving extra dimensions 4^ their masses and 
mixing may depend on energy through these running effects. Usually the scale of new physics 
in extra dimensional models is around or above the TeV energy scale. In order to be relevant 
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for the LSND/MiniBooNE problem one has to assume that the mechanism responsible for 
the power law running can be extended to the energy scale relevant for the experiments 
under consideration (MeV to GeV) in the active-sterile neutrino sector. 

At this point I will not go into further details and leave the question whether indeed a full 
model for MED or EDM oscillations can be constructed from any of the mentioned mecha- 
nisms open for future work. I add that in a given realisation the energy dependence might 
be different than assumed in Eq. ([T]). However, the generic assumption of a power law should 
be a reasonable approximation in many cases and capture the relevant phenomenology. 

VI. CONCLUSIONS 

I have considered short-baseline (SBL) neutrino oscillation data including LSND and 
MiniBooNE in the framework of sterile neutrino oscillations, assuming that the properties 
of the sterile neutrino depend on its energy in a rather exotic way. Along these lines I 
considered two different scenarios. First, I have assumed that the mass of the sterile neutrino 
scales with its energy as (0 < r < 1). This introduces a modified energy dependence 

(MED) in oscillations: Instead of the standard l/E^ dependence one obtains a MED with 
Second, I have assumed an energy dependent mixing (EDM) of the sterile neutrino 
with the active ones, namely that the elements of the mixing matrix |f/e4p and |f/^4p scale 
with 1/El,. In a given model realisation one can expect that both, masses as well as mixing 
depend on energy in a correlated way. Here I have not specified any underlying theory, and 
a phenomenological analysis has been performed assuming the presence of either MED or 
EDM separately, to show the impact on the global fit of SBL data. For a given model it is 
easy to generalise the analysis and estimate the effect of the simultaneous scaling of masses 
and mixing. 

I find that under the hypothesis of MED or EDM oscillations LSND and MiniBooNE 
data become fully consistent, and the bound from disappearance data overlaps with the 
LSND allowed region. The global fit including all relevant appearance and disappearance 
experiments improves by 12.7 (MED) or 14.1 (EDM) units in with respect to the standard 
(3+1) oscillation case, and the consistency of LSND with no-evidence appearance experi- 
ments and with disappearance experiments improves from 4 x 10~^ for standard oscillations 
to 3% (MED) or 5% (EDM). If the tension between LSND and KARMEN is removed from 
the analysis perfect consistency of all data is found, with probablities of 18% and 35% 
for MED and EDM, respectively, whereas the probablitiy of standard oscillations remains 
unacceptably low. Consistency of the global data is obtained in the MED framework by 
shifting the sensitivity of high energy experiments like MiniBooNE and CDHS to larger 
values oi L/ with respect to low energy experiments like LSND and Bugey. In the case of 
EDM oscillations the sensitivity of high energy experiments to the mixing angle gets weaker 
compared to low energy experiments, leading to consistency of all data. 

In summary, I have shown that under the assumption of a non-standard energy de- 
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pendence of sterile neutrino oscillations the description of global SBL data is significantly 
improved. This result is based on the fact that various experiments operate at different 
energy regimes. 
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